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SPECIFIC  HEAT  OF  VITREOUS  SILICA 


Abstract 

A nethod  is  developed  for  the  study  of  the  thermal 
motions  of  the  atoms  in  vitreous  silica.  This  enables  ten- 
tative calculations  of  specific  heat  based  on  various  forms 
of  binding  between  atoms  and  on  vibrations  of  single  atoms 
independent  of  their  neighbors.  First  comparisons  with  experi- 
ment are  not  very  satisfactory  and  methods  are  developed  for 
the  use  of  larger  and  larger  groups  in  order  to  approach  more 
closely  to  the  actual  network  vibrations. 

Introduction 

In  an  earlier  paper^^it  was  shown  that  a specific  heat 
curve  agreeing  reasonably  well  with  the  experimentally  deter- 
mined curve  cculd  be  calculated  from  a simple  model.  This 
model  consisted  of  silicon  atoms  which  could  oscillate  in 
three  directions  with  three  equal  frequencies  and  oxygen 
atoms  which  could  oscillate  in  two  transverse  directions  be- 
tv:een  pairs  of  silicon  atoms  with  two  equal  frequencies  and 
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with  a different  higher  frequency  longitudinally  between  the 
silicon  atoms . The  values  of  the  frequencies  were  not  derived 
from  any  detailed  model  of  the  structui’e  of  vitreous  silica 
but  were  rather  chosen  to  give  the  best  fit  with  the  experi- 
mental curve.  It  seemed  worth  while  to  extend  this  work  by 
trying  to  derive  the  curve  of  specific  heat  versus  temperature 
for  vitreous  silica  from  our  knowledge  of  its  structure  and 
some  reasonable  assumptions  as  to  the  nature  of  the  forces 
between  the  various  atoms  making  up  the  structure. 

Nature  of  Interatomic  Forces 

Considerable  doubt  exists  as  to  what  is  a reasonable 
assumption  as  to  the  nature  of  these  forces.  On  the  one  ex- 
treme there  might  be  considered  the  purely  ionic  picture 
where  the  forces  would  be  Coulomb  attractions  or  repulsions 
modified  only  by  repulsions  varying  with  a power  of  the  se- 
paration much  higher  than  the  second.  In  this  case,  in  com- 
puting the  field  of  force  in  which  each  ion  moves  one  must 
consider  the  effects  arising  from  a large  number  of  neighbors 
since  the  Coulomb  forces  fall  off  so  slowly  with  distance. 

On  the  other  extreme  the  forces  might  be  those  arising 
from  pure  covalent  binding.  In  this  case  the  interatomic 
interactions  would  arise  only  from  close  neighbors  but  some 
doubt  would  arise  as  to  the  exact  mathematical  form  to  use 
to  express  these  interactions.  For  purely  central  forces  a 


( 


(2) 

Morse  function  night  bo  choeen  such  as  to  give  the  correct 
heat  of  fomation  but  the  effect  of  specially  directed  covalent 
binding  might  still  be  difficult  to  handle. 

In  all  probability  the  true  picture  is  somewhere  between 
these  two  extremes,  Pauling,  'on  the  basis  of  the  electo- 
negativities  of  silicon  and  oxygon,  describes  the  binding  as 
50#  ionic.  In  order  to  obtain  more  experimental  indication 
as  to  the  nature  of  the  binding  it  is  justifiable  to  consider 
either  of  the  two  extremes  and  use  the  lack  of  agreement,  if 
found,  as  a stepping  stone  to  a more  exact  description.  The 
present  paper  considers  mostly  the  ionic  model. 

Ionic  Model  of  Vitreous  Silica 

As  in  a previous  paper  by  Smyth,  Londeroe  and  Lorey,^ 
each  silicon  ion  is  considered  to  be  contained  in  the  field 
arising  from  four  oxygens  tetrahedrally  arranged  around  the 
silicon,  four  silicons  at  twice  this  Si-0  distance  and  lying 
along  the  extension  of  the  Si-0  lines  and  twelve  oxygens  at 
a slightly  greater  distance.  Only  one  half  of  the  effect  of 
these  twelve  oxygens  is  considered  so  that  tho  charge  on  those 
ions  affecting  each  silicon  ion  is  -4  and  the  remaining  structure 
outside  these  is  considered  as  neutral  and,  because  of  the 
increasing  randomness  of  the  glass  structure  at  these  separations 
giving  effects,  which,  averaged  ovor  tho  whole  structure  will 
be  effectively  zero. 
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Similarly  each  oxygen  ion  i3  contained  in  the  field  of: 

(1)  Two  Si  ions  diamotrically  opposed  to  each  other. 

(2)  Six  oxygens.  (The  remaining  nearest  neighbors  of 
the  above  two  silicon  ions). 

(3)  One  fourth  of  the  six  silicon  ions  at  the  second 
closest  Si-0  distance. 

These  surroundings  exactly  balance  the  charge  of  -2  on 
each  oxygen  ion.  Tho  effect  of  tho  remainder  of  the  structure 
is  neglected,  being  neutral. 

The  question  of  the  validity  of  neglecting  these  neutral 

surroundings  is  one  which  deserves  a more  rigorous  treatment 

beyond  the  scope  of  this  paper.  In  crystals,  because  of  tho 

regularity  of  structure  it  cannot  be  neglected  but  must  be 

(5') 

exactly  considered,  as  in  the  methods  of  Dwald  and  others. 
The  assumption  is  made  here,  however,  that  the  lack  of  a re- 
gularly repeating  structure  beyond  the  neighbors  considered 
justifies  our  neglect  of  any  contribution  to  the  electrostatic 
field  from  thi3  source.  Any  conclusions  reached  must  there- 
fore be  subject  to  critical  examination  from  the  standpoint 
of  this  particular  assumption. 

Calculation  of  Frequencies 

As  a first  approach  to  tho  problem  of  studying  tho 
vibrations  of  the  atoms  it  might  be  considered,  as  in  tho 
Einstoin  theory  of  specific  heats  of  crystals,  that  each  atom 


vibrates  in  a manner  independent  of  the  motions  of  its  neigh- 
bors and  in  the  fiold  which  would  arise  from  all  of  its 
neighbors  being  in  their  equilibrium  positions.  Later  there 
shall  be  taken  up  the  problem  of  allowing  also  for  the  effects 
of  the  vibrations  of  each  atom  on  its  neighbors. 

Tho  restoring  force  acting  on  any  atom  when  it  is  dis- 
placed a distance  in  any  direction  from  its  equilibrium 
position  is  the  same  in  direction  and  magnitude  as  if  it  had 
been  at  rost  and  all  of  the  other  atoms  had  been  moved  a distance 
■ii  in  the  opposite  direction.  The  first  problon  is  therefore 
to  compute  the  force  acting  on  an  ion  at  the  origin  if  an  ion 
of  known  charge  at  somo  location  (a,  b,  c)  is  displaced  to  a 
new  nearby  position  (a+u,  b+v,  c+w).  This  must  tako  into 
account  both  tho  Coulomb  forces  and  also  the  higher  power 
repulsions  which  are  important  for  very  near  neighbors. 

As  far  as  the  Coulomb  forces  are  concerned  this  is 
equivalent  to  finding  the  eloctric  field  at  the  origin  arising 
from  a dipole  placed  at  (a,  b,  c)  whose  axis  is  parallel  to 
the  displacement  (u,  v,  w)  and  whoso  magnitude  is  the  charge 
on  the  ion  at  (a,  b,  c)  multiplied  by  tho  linear  displacement. 
This  is  a well  known  expression  worked  out  in  any  text  book 
on  electromagnetic  theory. ^ Putting  it  in  tho  form  con- 
venient to  tho  problem  in  hand  tho  expression  may  be  writton: 
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where  Fx,  Fy  and  Fz  aro  respectively  the  x,  y and  z components 
of  the  force  on  an  ion  of  charge  Oj.  at  the  origin  due  to  dis- 
placement (u,  v>  w)  on  an  ion  of  charge  ca  located  at  (a,  b,  c), 
and  */"  is  tho  distance  from  the  origin  to  the  pont  (a,  b,  c)  so 
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It  might  bo  montioned  that  those  relations  hold  rigidly 
only  when  tho  displacements  arc  small  comparod  with  the  intcrionic 
distances  involved. 

The  forces  arising  from  tho  highor  power  repulsions  can 
bo  expressed  in  very  much  tho  samo  form.  If  tho  ion  at  tho 
origin  and  the  ion  at  (a,  b,  c)  ropcl  each  other  with  a forco 
of  magnitude 
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whon  separated  a distance  '/*  from  each  other  then  the  components 
of  tho  forco  on  the  ion  at  the  origin  when  the  ion  et  (a,  b,  c) 
is  displaced  an  amount  (u,  v,  w)  aro: 
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Here  tho  approximation  is  not  quite  so  good  and  tho 
neglecting  of  higher  powers  of  u,  v and  w than  tho  first  be- 
comes less  and  loss  justifialbc  as  Kh  increases. 

The  modol  to  bo  used  in  calculating  tho  natural  froquoncies 
of  each  ion  can  now  be  sot  up  and  the  mathematical  exprossions 
worked  out  for  those  frequencies. 

Ionic  Modol 


Tho  interactions  considered  in  this  treatment  will  be: 


(1)  Coulomb  Interactions  botwcon  oach  ion  and  those 


noighbors  previously  discus  sod. 

(2)  Ropulsions  of  tho  form  — botwoen  oach  silicon 


ion  and  its  closost  oxygen  noigjhbors. 

(3)  Repulsions  of  tho  form  —fir  botwoen  oach  oxygon 
and  its  closost  oxygon  noighbors, 

Ropulsions,  other  than  oloctrostatic  repulsions  will  not 
bo  considered  botwoen  silicon  ions  and  thoir  nearost  silicon 
noighbors  or  between  any  pairs  of  ions  with  groator  separations. 

In  ordor  to  computo  tho  froquoncios  of  a silicon  ion  tho 
following  arrangement  i3  considered.  Tho  silicon  ion  placod 
at  the  origin.  Its  neighbors  aro  locatod  as  follows: 

(1)  Four  oxygon  ions  at  (s,  s s),  (s,  -s,  -s),  (-s,  s,  -s) 
and  (-s,  -s,  s),  Theso  aro  totrahedrally  arranged  about  tho 
silicon  ion  at  a distance  such  that 

3b1 

(2)  Four  silicon  ions  at  (2s,  2s,  2s),  (2s,  -2s,  -2s) 

(-2s,  2s,  -2s)  and  (-2s,  -23,  2s), 

Those  arc  at  a distanco  t such  that  ~fc,  * <-  , 

‘•e  ' 

(3)  Twolve  oxygen  ion3  at 
(3s,  3a,  a)  (3s,  a,  3s)  (s,  3s,  33) 

(-3e,  -3s,  s)  (-3s,  -s,  3s)  (-s,  -3s,  3a) 

(3s,  -3s,  -s)  (3s,  -s,  -3s)  (s,  -3s,  -3s) 

(-3s,  3s,  -s)  (-3s,  8,  -3s)  ( — s , 3s,  -3s) 

Tho  total  offect  of  theso  twolvo  will  be  computod  and  then 
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divided  by  two  ao  that  the  neighbors  altogether  will  just 
exactly  balance  the  charge  of  the  silicon  ion  at  the  origin, 
while  we  at  the  same  time  avoid  the  necessity  of  making  in- 
vidious distinctions  between  these  twelve  oxygen  ions.  Each 
group  of  three  of  these  oxygens,  together  with  one  of  the 
oxygens  of  group  (1)  tetrahedrally  surround  ono  of  the  silicons 
of  group  (2).  Each  of  these  oxygens  is  a distance  from  the 
origin  where 


All  of  these  neighbors  (1),  (2)  and  (3)  will  be  displaced 
an  amount  (u,  v,  w)  and  the  resultant  force  on  the  silicon 
ion  at  the  origin  computed. 

Using  equations  1A,  IB  and  1C  it  can  be  easily  shown  that 
if  all  of  the  oxygens  at  distance"^  from  the  origin  are  given 
the  same  small  displacement  (u,  v,  w)  the  electostatic  forces 
alone  give  a zero  resulting  force  on  the  silicon  ion  at  the 
origin.  Also  if  all  of  the  silicons  at  distance  V^cr  all  of 
the  oxygons  at  distance  arc  disolacod  by  the  same  displace- 
ment (u,  v,  w)  the  resulting  force  on  the  silicon  ion  at  the 
origin  is  zero.  Hence,  if  all  of  these  neighbors  are  held 
at  rest  and  the  silicon  ion  is  given  a small  displacement 
there  is  no  restoring  force  arising  from  electrostatic 
attractions  or  repulsions.  All  of  the  restoring  force  on 
the  silicon  ion  arises  from  the  higher  power  repulsive  forces 
of  its  neighboring  oxygens  (if  we  can  neglect  the  effect  of 
the  high  power  repulsions  of  more  distant  ions). 
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Using  equations  i|A,  i|B  and  4C  it  can  then  bo  shown  that 
if  the  four  oxygens  in  group  (1)  are  given  each  a displacement 
of  (u,  v,  w)  then  the  resulting  force  components  on  a silicon 
at  the  origin  will  be 
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from  which  it  can  immediately  be  deduced  that  if  a silicon 
ion  is  given  a displacement  AL  in  any  direction  it  experiences 
a restoring  force 

a ■ > h , 

-Jr 

in  the  exactly  opposite  direction.  The  field  in  which  it 
moves  is  quite  spherically  symmetrical  for  small  displacements. 
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Forces  on  Oxygen  Ions 

In  the  case  of  0x73 ur.  ions  more  attention  has  to  be  oaid 
to  the  direction  of  displacement.  It  i3  obvious  that  there 
is  likoly  to  be  a difference  in  the  restoring  force  depending 
on  whether  the  displacement  is  along  the  line  of  the  two 
neighboring  silicon  ions  or  transverso  to  it.  It  will  bo 
shown  that  if  throe  mutually  perpendicular  directions  aro 
chosen  such  that  one  of  them  lies  along  the  silicon-silicon 
direction  then  displacements  and  vibrations  in  any  one  of 
these  directions  are  independent  of  what  is  happening  in  the 
other  two.  This  docs  not  hold  for  three  mutually  perpendicular 
directions  not  so  chosen.  The  model  will  be  3ct  up  as  follows: 

The  oxygen  ion  in  question  will  be  located  at  the  origin. 
Its  neighbors  will  be: 

(1)  Two  silicons  at  (a,  s,  s)  and  (-3,  -3,  -a) 

(2)  Six  oxygens  at  (2s,  2s,  0)  (2s,  0,  2s)  (0,  2s,  2s) 
(-2s,  -2s,  0)  (-2s,  0,  -2s)  (0,  -2s,  -2s),  These  are  at  a 

r — 

distance”^  where  ^ ^ or  ~ \ *5 

(3)  Six  silicons  at  (3s,  3s,  -s)  (33,  -s,  3s)  (-s,  3s,  3s) 
(~3s,  “3s,  s)  (*3s,  s,  -3s)  (s,  -3s, “3s).  Of  this  last  group 
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we  divide  the  total  effect  by  four  for  the  sane  reasons  discussed 
previously. 

In  order  to  study  longitudinal  vibrations  (i.e.  along  the 
silicon-silicon  line)  a displacement  of  (u,  u,  u)  will  be  im- 
posed in  place  of  the  general  displacement  (u,  v,  w)  used  for 
the  silicon  ion. 

A displacement  of  (u,  u,  u)  imposed  on  the  two  silicon 
ions  in  group  (1)  gives  rise  (using  equations  1A,  IB  and  1C) 
to  an  electrostatic  force  on  the  oxygen  at  the  origin  whose 
components  are: 
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Therefore  a displacement  of '-CL  along  the  silicon-silicon  line 
of  the  oxygen  ion  at  the  origin  brings  to  bear  on  it  an  electro- 
static force  of 


~~hi 


(6) 


-13- 


meaning  that  as  far  as  the  electrostatic  attractions  of  the 
two  silicons  alone  are  concerned  the  oxygen  would  be  in  un- 
stable equilibrium. 

Using  equations  (4A,  4B#  and  I4.C ) the  force  on  the  oxygen 
at  the  origin  arising  from  higher  power  repulsions  of  the  two 
silicons  is 
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for  e displacement  It  of  this  oxygen  along  the  silicon-silicon 
line.  This  leads  to  stable  equilibrium. 

A displacement  “i/.  of  the  oxygen  ion  at  the  origin  along 
the  silicon-silicon  line  will,  because  of  eloc trostatic  forcos 
from  the  feix  oxygens  in  group  (2)  give  rise  to  a force  on  tho 
oxygen  ion  of 
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The  repulsive  forces  of  these  six  oxygens  give  a force  on  the 
oxygen  moved  longitudinally  a distance^  of 
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Tho  six  silicons  in  group  (3)  wi] 1 cause  a force  (purely  electro- 
static) on  the  oxygen  when  moved  a dis tance v / /.  along  the 
silicon-silicon  line  of 
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Again  the  two  groups  of  attracting  ions  pulling  against 
each  other  tend  to  produce  unstable  equilibrium.  Of  this 
quantity  only  one  quarter  is  considered  for  reasons  already 
discussed  leaving 
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Using  the  relations  that 
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tho  complete  effect  may  be  summed  tip.  The  resultant  force 
on  an  oxygen  ion  displaced  a distance yUL  along  the  silicon- 
silicon  line  is 
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or  arranging  it  all  in  terms  of  ~Lt  the  nearest  silicon  oxygen 
distance,  and  going  through  the  various  numerical  calculations 
the  final  result  is 
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It  so  happens  that  the  second  and  third  torn  together  are 
numerically  greater  than  the  first  resulting  in  stable  equili- 
brium. The  restoring  force  increases  as  is  decreased. 

Transverse  Oxygen  Vibrations 

To  study  vibrations  of  the  oxygon  ion  in  a direction 
perpendicular  to  the  silicon-silicon  line  each  other  ion  i3 

given  a displacement  perpendicular  to  this  line.  The  direction 

> / ' 

cosines  of  the  silicon-silicon  line  ero  (X , ^ ) . A 

direction  with  diroctions  cosinos  (1,  m,  n)  will  bo  perpendicular 
to  this  if  1 + in  + n B 0.  If  therefore  each  ion  is  given  a 


displacement  l 


— | where  u and  v are  arbitrary  this 

displacement  will  always  be  perpendicular  to  the  longitudinal 
direction  and  proper  choice  of  the  ratio  of  -U.  to  IT  will  cover 
all  possible  such  transverse  diroctions.  It  will  be  shown  that 
in  each  case  the  restoring  force  will  make  an  anglo  of  0 or  l80° 
with  the  displacomont . This  holds  only  for  .displacements  parallel 
to  or  perpendicular  to  the  longitudinal  direction  and  would 
not  hold  for  any  direction  in  between  (except  for  the  special 
case  whothor  tho  rostoring  forces  in  the  longitudinal  and 
transverse  diroctions  are  the  samo). 

First  tho  two  silicons  at  (s,  s,  s)  and  (-s,  -s,  -s)  are 
given  a displacement^  , From  equations  4A,  Lj.B 

and  ]+C,  the  components  of  the  forco  on  tho  oxygen  ion  ett 
(0,  0,  0)  aro 
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Therefore  if  the  oxygen  ion  is  given  a displacement m a 
transverse  direction  there  is  exerted  on  it  a force 
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which  means  a positive  restoring  force  back  along  this  same 
direction. 

The  repulsions  of  these  two  silicons  give  rise  to  a force 


(17) 


if  the  oxygen  is  given  a displacement '~iL  in  a transverse  direction. 

Similarly  the  six  oxygens  of  group  (2)  by  electrostatic 
action  give  a force  of 
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on  an  oxygon  ion  at  the  origin  moved  transversely  a distance  -CL, 
The  non-Coulonb  repulsions  of  these  oxygens  give  a force  of 
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Taking  ono  quarter  of  the  effect  of  the  oloctrostatic 
offset  of  the  six  silicons  at  there  is  obtained  for  tho 
forco  acting  on  anoxygen  at  the  origin  displaced  transversely 
a distance  LL 
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Tho  complcto  transverse  effect  may  then  be  summed  up  to 


give 

lit 


iZt 

h 


-zf*  'it  + - 


•» Hi  + &*+> ({f-iCyit 

1 ri  v.  / 


Putting  this  in  terms  of  only  there  is  obtained 


\\  f i 


~z  h 


yp • < p -»vj -it 


-t: 


Stability  depends  on  tho  first  torn  boing  numorically 
greater  than  tho  sum  of  tho  othor  two.  The  restoring  force 
here  decreases  as  docreasos. 
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Choloe  of  Constanta 

There  remains  the  problem  of  selecting  the  most  probable 

( 7) 

values  for  A,  B and  W According  to  Fowler'  'two  neon  like 
ions  might  be  expected  to  repel  each  other  with  a force  varying 
as  the  inverse  eleventh  power  of  the  separation  and  v'uls  there- 
fore taken  as  eleven.  Furthermore  there  is  listed  also  in 
Fowler  a value  for  B giving  the  magnitude  of  the  interaction 
between  two  oxygen  ions.  This  value  is  if  -f  is  measured  in 
centimeters  and  the  force  in  dynes. 


3 - /?•  6 * ID 


(23) 


A must  then  be  chosen  so  that  the  energy  of  the  network  is 
a minimum  for  the  known  value  of  the  silicon-oxygen  separation 
in  vitreous  silica. 

Energy  of  Network 

Limiting  attention  to  interactions  between  each  ion  and 
those  neighbors  previously  considered  it  may  be  easily  shown^ 
that  the  electrostatic  contribution  to  the  energy  of  the  network 
is 

_ 3^  t!"  ' 3 2 £ L 

r -A 


(24) 
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which  is  a quantity  becoming  increasingly  negative  as  "h 
(along  with  , 'j^and  ”)4  in  the  same  ratio)  is  decreased. 
This  reduces  to 

- /') . 6 I 4- 

1 I 

This  negative  energy  is  offset  by  the  positive  energy 
arising  from  the  repulsions  between  ions  as  they  come  close 
together.  For  the  structure  considered  this  is 


A.  m 
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v. 
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C.  • - I 
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The  total  energy  is  therefore 


(:& 


>v"' 


(27) 


The  equilibrium 

dis  tance 

' } , will  be 

such  that  this 

energy 

is  a minimum  or 

that  its 

derivative 

with  respect  to 

-'L, 

will  be  zero. 

Hence 

X. 

4- Pi 

n “ 1 

/'j\  i 1 - 
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0 from  which 


i 


,i  --  -3.  '102  cV 
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vs 
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/7<; 

1 i.  63 

(29) 

|i 

1 

i 

1 

Us  ing 

1 

< 

c 

= 4.77  x 

10"10  l.  b 'IL 

a 

S 

a 

1 

t; 

= 1.62  x 

10“8 

i 

\ 

i 

3 

= 17.6  x 

10-68 

3 

* 

= 11 

3.90200 

- .682  x 

10-88 

i.$  B 

= .196  X 

10-88 

■33 

Hence 

/) 

= .436  X 

10-88 

(30) 

Calculation  of  Frequencies 

Equations  (6),  (15)  and  (22)  may  now  be  used  to  compute 
in  actual  numerical  values  the  restoring  forces  on  the  oxygen 
and  silicon  ions,  and  hence  calculate  the  natural  frequencies 


of  oscillation  of  each  of  the  three  different  modes  of  vibration. 


(. 
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For  a silicon  ion  displaced  in  any  direction  the  restoring 
force  is 


178.50  x 104  dynes  per  cm.  (31) 

For  a longitudinal  displacement  of  an  oxygen  ion  tho 
restoring  force  is 


247.13  x 104dynes  per  cm. 


(32) 


For  a transverse  displacement  of  an  oxygen  ion  it  is 


50.75  x 104  dynes  per  cm. 


(33) 


If  the  restoring  force  per  unit  displacement  of  a 
mass  hi  is  then  the  frequency  of  oscillation  is 


V 


J ! -/>, 

b V V> 


(34) 


Using  the  knovm  masses  of  the  silicon  and  oxygen  ions, 
the  frequency  is  easily  computed.  The  relation  connecting 
frequency  with  the  characteristic  temperature  ,0  to  be  used 
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in  the  Einstein  specific  heat  functions  is 


(35) 


where  h is  Planck's  constant  and 'li.  is  Boltzmannte  constant. 
Table  I lists  these  numbers  for  the  three  modes  of 
vibration. 


Table  I 


Ion 

Silicon 

Oxygen 

Oxygen 


Direction  of  Vibration 
All 

Longitudinal 

Transverse 


Characteristic 

Temperature 

6L  1495° 
fcL  2328° 
t\  1054° 


The  specific  heat  can  be  calculated  from  these  characteristic 
temperatures  using  the  relation 


where  £ 


is  the  usual  Einstein  function 


ZroO 


'L  L 

l £ 


(^"0 


(37) 
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R is  the  gas  constant  and  M is  the  molecular  weight  of  silica. 

O 

The  specific  heat  computed  in  this  way  is  shown  plotted  in  fig.  1° 
along  with  the  experimental  values.  The  agreement  is  not  good, 
and  the  theory  must  be  examined  to  find  where  it  should  be 
modified  to  give  better  agreement  with  experiment.  By  comparison 
with  values  of  the  characteristic  temperatures  which  have  to 
be  assumed  to  get  reasonably  good  agreement  the  characteristic 
temperatures  used  here  are  all  too  high  which  means  the  fre- 
quencies also  are  too  high. 

Before  raising  the  question  of  the  validity  of  the 
completely  ionic  model,  it  is  expedient  to  examine  how  the 
Einstein  treatment  of  the  specific  heat  can  be  modified 
to  permit  an  approach  toward  the  Debye  method  while  retaining 
at  the  same  time  the  connection  between  the  calculations  and 
the  atomic  model.  The  characteristic  frequencies  have  been 
calculated  on  the  assumption  that  each  ion  vibrates  in  a 
manner  independent  of  the  vibrations  of  its  neighbors  which 
is  obviously  not  the  case.  If  an  ion  is  displaced  from  its 
equilibrium  the  effect  at  each  other  ion  is  as  if,  to  a first 
approximation,  there  were  added,  at  the  location  of  the  dis- 
placed ion  a dipole  of  moment  equal  to  the  charge  on  the  ion 
multiplied  by  its  displacement.  Each  vibrating  ion  therefore 
creates  a harmonically  varying  field  at  the  location  of  each 
other  ion  of  a magnitude  and  direction  which  depends  on  the 
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geometry  of  each  pair  of  ions.  In  very  much  the  same  way  as 
equations  1A,  IB  and  1C  were  derived  it  can  be  shown  that  if 
an  ion  of  charge  ee  at  the  point  (a,  b,  c)  is  displaced  a 
distance in  a direction  whose  direction  cosines  are 
(la,  nva,  rie  ) then  there  is  exerted  on  an  ion  of  charge  e^ 
at  the  origin  a force  whose  component  in  the  direction  (lx,  mx , 
ih  ) is: 


'Ci it 


where 


'L  . 0_ 

Cl  Xr  -r  (L> 


The  effect  of  the  repulsion  terms  can  be  considered  in  much 
the  same  way  and  the  result  is  that  if  two  ions  distance  ' 
from  each  other  repel  each  other  with  a force 


(39) 


then  the  component  of  force  along  the  (lx,  mx  , i>x  ) direction 
exerted  on  tho  ion  at  the  origin  by  a displacement  ^ of  tho 
ion  at  (a,  b,  c)  in  a direction  (la,  ith  , nx  ) is 


(40) 


Instead  then  of  handling  one  ion  independent  of  the  motions 
of  its  neighbors  one  may  then  select  a group  such  as,  for 
instance,  a single  SiOa  group  or  two  SiOE  groups  tied  together 
or  a group  of  say,  five  SiOa’s  arranged  in  a compact  unit 
and  study  the  vibrations  of  one  of  these  group  aa  a whole 
instead  of  merely  oach  isolated  ion.  In  each  case  the  ions 
outside  the  group  under  study  will  be  considered  as  being  at 
rest.  The  frequencies  will  then  bo  the  frequencies  of  the 
independent  modes  of  vibration  of  the  group  as  a whole  and 
in  this  way  one  might  expect  to  find  some  modes  with  con- 
siderably lower  frequencies  than  those  of  tho  ions  vibrating 
in  fixed  surroundings  which  would  allow  the  specific  heat 
curve  to  rise  more  rapidly  at  low  temperatures.  If  the  group 
were  made  large  enough  and  the  assumptions  as  to  interatomic 
forces  were  correct  then  one  would  expect  to  be  able  to 
calculate  correctly  the  specific  heat.  Tho  hope  is  that, 
because  of  the  lack  of  similarity  of  neighbors  beyond  the 
fourth  or  fifth  nearest  neighbors  the  groups  chosen  could  be 
small  enough  for  reasonable  ease  of  calculation.  If  this 
should  prove  to  be  the  case  we  would  have  a useful  tool  for 
oaloulating  other  physical  properties  such  as  coefficient 
of  expansion.  Free  energy  at  any  tompoisture  could  be  computed 
by  the  regular  methods  of  statistical  mechanics  and  that  net- 
work spacing  chosen  at  oach  temperature  such  that  the  free 


energy  is  a minimum.  One  migfrt  even  be  optimistic  enough 
to  foresee  that  the  effect  of  additions  such  as  any  of  the 
usual  network  modifiers  might  be  so  computed. 

As  an  illustration  one  can  take  the  case  of  a single 
SiOa  group,  i.o.  any  silicon  ion  and  any  two  of  it3  immediate 
oxygon  neighbors.  There  will  be  nine  degrees  of  freedom 
and  thero  should  therefore  be  nine  modes  of  vibration  of  this 
group.  The  Si  will  bo  taken  at  (0,  0,  0)  one  oxygon  at 
(s,  s,  8)  and  the  other  at  (-s,  -s,  s).  Any  three  mutually 
perpendicular  axes  may  be  choson  for  each  ion  and  one  must 
then  investigate  what  linear  combinations  of  those  will  give 
those  nine  normal  co-ordinates  which  will  represent  the  nine 
independent  modes  of  vibration. 

In  order  that  the  three  directions  chosen  for  an  oxygen 
ion  be  independent  of  one  anothor  one  should  be  along  the 
oxygen-silicon  line  and  the  other  two  transverse  to  it.  It  la 
also  convenient  to  choose  one  of  these  transverse  directions 
for  each  oxygen  perpendicular  to  the  plane  of  the  S10a  group. 
For  the  silicon  the  throe  directions  of  motion  will  be  chosen 
(1)  porpondicular  to  the  piano  of  the  SiOa  group,  (2)  along 
the  bisector  of  the  O-Si-O  angle  and,  (3)  in  the  Si08  plane 
perpendicular  to  tho  bisector.  The  direction  cosines  of  these 
directions  are  given  below: 


First  oxygon: 
Co-ordinato 


Q 

Co-ordinate  xa 

Co-ordinate  Xs 
Second  oxygen 
Co-ordinato  x« 

Co-ordinate  x^ 

Co-ordinate  x^ 

Silicon 
Co-ordinate  Xy 

Co-ordinate  xg 

Co-ordinate  x^ 

If  co-ordinate  xm  is 
exerted  on  the  ion  whose 
diroctioh  is  Fpjp  then  the 
that 

i— 

c 


-27- 


"28- 


Prom  the  symmetrical  form  of  equations  (38)  and  (l+O) 
it  is  apparent  that 

The  three  quantities  kgi  and  kL  and  k-^  will  be  defined 
as  follows. 

If  a silicon  ion  is  displaced  a distance  ^4  while  all 
other  ions  are  held  at  rest  the  force  on  the  silicon  ion  will  be 

ksitU. 

If  an  oxygon  ion  is  displaced  a distance  '-LC  along  the  oxygen- 
silicon  lino  the  force  exerted  on  it  will  be  (all  other  ions 
remaining  in  positions) 


*1  % 

If  an  oxygen  ion  is  displaced  a distance  ^ in  a transverse 
direction  the  force  on  it  will  be  (all  other  ions  remaining 
positive ) 


kT  -a 

For  equilibrium 
negative  quantities, 
given  in  expressions 


it  is  clear  that  k^^and  k^.  are  all 
Their  numerical  values  have  already  been 
(31),  (32)  and  (33). 
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Si  nee  a number  of  the  k^^s  go  to  zero  it  in  convenient 
to  evaluate  these  before  writing  down  the  equations  of  motion 
in  each  of  the  nine  co-ordinates.  Prom  our  choice  of  co- 
ordinates it  is  clear  that  k12,  k23 , k31,  k^,  k^6,  k^,  k?g, 
kg^,  and  kgy  are  all  zero.  The  values  of  the  rest  are  given 
below; 


I ii  n-  *■  i ) 


0 
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The  nine  equations  of  motion  may  then  be  set  up  omitting 


terms  for  which  the  k*s  are  zero. 


KjL  Lf  + kst4  ^4.  4c  -A 

A,  vi^ 

-Wc  2F- 

(44) 

“4  ** 

*7  - ;UC  - 

(45) 
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U L 3? 
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(46) 

i+A.K  = 

rCv 

»l c TP- 

(47) 

P + 

5 k-3  9 4 

m ^ 
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4 *->  * 

- '4-  £t  f 

4>v*U 

Cx.t 
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(49) 

A,  ^ + 

4; 

i«+  Ol 

- 11  j« 

d\7 

>c  ^?!r  i 

(50) 

/6/s.  i'.,  +■ 

A*  *-3 

Cl, 
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(5D 
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*-%,k+ 

K\+As 
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(52) 

Co-ordinates  Xa , xg  and  xy  are  independent  of  all  the 
other  co-ordinates.  Further,  because  k^y  = k^y  equations 
(45)  and  (46)  can  be  combined  to  give 
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c 


t- 


) 


(S3) 


( 


-33- 


which  gives  immediately  one  of  the  independent  modes  of 
vibration  and  corresponds  to  both  oxygen  ions  vibrating 
perpendicular  to  the  SiOa  plane  so  that  one  is  always  as 
far  on  one  side  as  the  other  is  on  the  other  side.  Such 
a motion  can  continue  without  exciting  any  of  the  other  co- 
ordinates into  motion. 

Similarly  equations  (45),  (48)  and  (50)  may  be  combined 
by  multiplying  both  sides  of  (50)  by  some  number  >4  and  adding 
(using  the  fact  that  k = k„  ) 

2?  57 


vr- 


(55) 


-17  -r  ^ '^Sc 

then  equation  (54)  will  reduce  to  a simple  harmonic  equation 
for  the  co-ordinates 


t2  + "yn  ^ 


(56) 


Actually  equation  (55)  yields  two  such  values  of  ^ giving 
all  of  the  three  (counting  equation  (53))  independent  modes 


C 
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1 

'♦ 

t 

I 


I 


-34- 


© 


of  the  co-ordinates  x$ , x and  x^. 

To  obtain  the  other  six  modes  of  vibration  it  should 


ons  ( 43 ) 

that 

kl6  = 

k34 

(57) 

kl8  = 

CO 

M* 

1 

(58) 

“l9  = 

k49 

(59) 

k36  “ 

"k68 

(60) 

k39  = 

k69 

(61) 

Using  these  relations  equations  (44),  (46),  (47),  (49) 
and  (52)  may  be  combined  to  give 


d*  \ 

('t-1  w (62) 

At(t*  r,)  f T6r(i, ->£,)+ V m° Jr  (*»* «*)  (63) 

^ ('t ♦ h)  f-  'K (u < t.)  6 4. 4 ‘ h\.  %j2  (64) 
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of  the  co-ordinetes  Xo , x^  and  x^« 

To  obtain  the  other  six  inodes  of 


be  noted  from  equations  (43)  that 


vibration  it  should 


kl6 

“ k34 

(57) 

PV 

H* 

CD 

"k48 

(56) 

^9 

k49 

(59) 

k38 

= ~k68 

(60) 

k39 

= k69 

(61) 

Using  these  relations  equations  (44),  (4&),  (47),  (49) 
and  (52)  may  be  combined  to  give 

q)  f 7$  ,-(4 ->£*)  + V 11 W-  (^  *■<■)  (63) 

(L+  l,  j f-  iy,  £*}  t-  ~k^  =■  ^ ^ ^r£  (61*) 
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If  (63)  is  nmltiplied  by-*/,  and  (6I4.)  by  /}/  there  is  obtained 


rh  * L + ^ * 3.  'Wft*  + *♦)  + + 1,  kst  * 3,  *><  h 

\ ( '* 

3,^)  t.  = M0  tfijii.'  t.)  fd(^  &.)j  + 5,  »»Si 


<ttv 
(65) 


This  transforms  to  a simple  harmonic  equation  in  the  variable 

r>7  s 

*i  + x*  + ^Ua  + x6^  +/i/  Jtq  where  0 = if  the  following 

two  relations  hold 


and 


k.L  -h  "kl4.  + *1  ^ ^ 
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These  may  be  rearranged  to  yield 
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t 


(69) 


are 


which  car.  be  solved  graphically  for  and  ,3 \ » There 
three  pairs  of  values  for  and  $ ^ yielding  three  more  of 
the  nine  modes  of  vibration.  The  final  three  may  be  found 
by  combining  (44),  (46),  (47),  (49)  and  (53 ) to  yield 

"kL('ir  A* 

--A*  (*,-«♦)*  'ir  ni) 


d\< 


(72) 


'A?  ^1"  ^4)  ^ ' ^Si.  c6£ 

Again  multiplying  (71)  by  0^  and  (72)  by jOx.  and  combining 
there  is  obtained  an  equation  which  will  be  a simple  harmonic 
equation  in  the  variable 

\)Lt' -h  06/  ^ £s~  0)L<z  - isjr 

/ 

are  connected  by  relations  similar  to  (68)  and  (69)  and 
which  actually  are 

'-4/i^X  - ')  t A ~ Al  ~ ' '&/*) 

'Il,^  /xL>  1 & 

>> 


3* 


■*■  $ o ^ 6,4 


-37- 


This  again  may  be  solved  graphically  yielding  three  pairs 
of  values  for  0u_  and  and  giving  the  final  three  modes  of 

i 

vibrations . 

It  remains  to  insert  the  numerical  values  for  the  k's 
and  determine  whether  the  move  toward  this  larger  grouo  is 
an  improvement  over  use  of  the  single  vibrations  of  the 
separate  ions.  Y/hen  this  is  done  there  are  obtained  nine 
frequencies  and  nine  characteristic  temperatures  for  the 
Si02  group.  These  temperatures  are  listed  below  in  order 
of  increasing  magnitude 

795°K 

962 

985 

1020 

1332 

1518 

1601 

2376 

2477 

Of  these  the  first  four  correspond  to  modes  where 
transverse  oxygen  vibrations  are  dominant.  The  next  throe 
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arise  from  modes  where  a silicon  vibration  is  dominant  and 
finally  in  the  last  two  a longitudinal  oxygen  vibration  is 
dominant.  The  larger  group  has  introduced  some  lower 
characteristic  temperatures  and  gives  a slightly  better  but 
still  not  good  agreement  with  the  experimental  values. 

Presumably  now  if  the  assumptions  as  to  structure  and  forces 
between  ions  we re  correct  and  if  the  various  approximations 
used  along  the  way  were  justified,  then  by  choosing  a large 
enough  neutral  group  of  ions  and  working  out  all  the  proper 
modes  of  vibrations  of  this  group  a specific  heat  curve  could 
be  computed  which  would  fit  the  experimental  data.  Unfortunately, 
to  work  out  mathematically  the  proper  modes  of  vibration  for 
any  group  larger  than  that  treated  above  requires  a prohibitive- 
ly long  computation. 

However,  it  is  possible  to  set  up  an  electrical  analogue 
in  which  the  movements  of  charges  arc  determined  by  exactly 
the  same  equations  of  motion  as  would  govern  the  movements 
of  ions.  As  an  illustration,  the  circuit  of  fig.  2A  re- 
presents the  circuit  for  an  ion  vibrating  in  one  direction 
independent  of  the  motions  of  its  neighbors.  The  charge  on 
the  condenser  corresponds  to  the  displacement  of  the  ion, 
the  inductance  L corresponds  to  the  mass  of  the  ion  and  the 
reciprocal  of  the  capacity  C corresponds  to  the  restoring  force 
per  unit  displacement. 
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By  virtuo  of  equation  (1+2)  this  analogue  may  be  extended 

to  two  ions  whore  tho  movement  of  one  causes  a force  on  the 

other,  and  vice  versa.  Here  represents  the  mas3  of  ion  (1) 

and  La  the  mass  of  ion  (2).  The  quantity  1,  + 1 represents 

Ci  Cc 

the  restoring  force  on  ion  (1)  because  of  all  its  neighbors 

at  rest  and  1+1,  represents  tho  corresponding  quantity  for 
Ca  Cc 

ion  (2).  1^  represents  the  force  on  ion  (1)  for  a unit  dis- 

Cc 

placement  of  ion  (2)  or  (by  equation  (42))  the  force  on  ion 
(2)  caused  by  a unit  displacement  of  ion  (1). 

Fig.  2c  represents  one  of  tho  possible  ways  in  which  throe 
degrees  of  freedom  could  be  coupled  togethor.  It  is  planned 
to  use  such  a system  in  tho  investigation  of  the  proper  modes 
of  vibration  of  larger  neutral  groups.  A large  froquoncy 
range  can  be  covered  with  a vacuum  tube  oscillator  and  the 
resonant  frequencies  can  easily  bo  picked  up  with  a vacuum 
tubo  voltmeter  or  oscillograph. 

Future  Work 

The  methods  developed  here  3ive  a means  of  quantitatively 
checking  the  experimental  values  of  specific  heat  with  any 
choson  field  of  f'orcc  between  the  constituent  atoms  cr  ions. 

The  first  agreement  is  not  good  and  the  most  obvious  places 
to  try  to  improve  the  assumptions  made  arc  as  listed  below: 
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(1)  Choice  of  puro  ionic  model.  Possibly  a mixture 
of  this  model  end  a model  showing  also  potential  energies 
between  nearest  neighbors  of  a fern  like  a Morse  function 
might  give  better  results.  Constants  would  be  chosen  so 

as  to  give  agreement  between  the  calculated  and  experimental 
heats  of  formation. 

(2)  Inclusion  of  more  ions  in  the  groups  considered  to 
vibrate  independently.  This  has  already  been  discussed. 

(3)  A good  detailed  geometrical  investigation  of  the 
probable  deviations  from  the  regular  cristobalitc  structure. 

It  should  be  possible  to  find  a limitod  number  of  ways  of 
introducing  the  irregularities  which  differentiate  the  vitreous 
state  from  the  crystalline  state.  From  the  results  of  the 
paper  on  Compressibility  of  Vitroou3  Silica  by  Smyth,  Londere 
and  Lorey  it  appears  that  bocause  of  such  irregularities 

some  but  not  all  of  the  Si-O-Sl  lines  are  actually  bent. 

This  would  lower  the  frequency  of  the  transverso  vibrations 
and  help  to  givo  a bettor  fit  with  tho  experimental  specific 
heat  curve  at  lower  temperatures. 

If  by  one  or  more  of  these  methods  it  is  found  possiblo 
to  get  a reasonably  good  agreement  between  calculated  and 
experimental  specific  heat  values  using  a reasonably  limited 
number  of  independent  modes  of  vibration  it  should  be  possible 
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to  sot  up  the  partition  function  for  tho  system  and  calculate 
tho  free  energy  at  any  tomporature  and  any  Si-0  spacing. 
Minimizing  tho  free  energy  at  oach  temperature  would  allow 
calculation  of  the  thermal  expansion.  In  view  of  the  fact 
that  at  moderately  low  tomperaturos  only  the  transverse  oxygon 
frequency  plays  any  important  part  and  sinco  it  is  easy  to 
show  from  equation  (22)  that  on  the  basis  of  tho  model 
considered  this  froquency  decreases  rather  than  increases 
with  decreasing  volume  it  seems  very  probable  that  this 
alone  will  be  sufficient  to  oxplain  the  nogative  and  low 
expansion  of  vitreous  silica  at  low  temperatures.  Extension 
of  thi3  work  to  glasses  containing  a second  oxide  might  not 
bo  out  of  tho  question. 

Conclusions 

The  results  show  that  a calculation  of  tho  frequencies 
of  vibration  of  silicon  and  oxygon  ions  separately  on  a pure 
ionic  model  is  not  sufficient  to  get  a good  agreement  between 
calculated  and  experimental  specific  heat  values  for  vitreous 
silica.  They  also  show  several  ways  in  which  this  very 
naivo  picture  may  be  modifiod  to  give  hotter  agreement. 

Tho  fact  that  tho  transverso  vibration,  which  is  the 
only  important  one  at  low  temperatures,  decreases  with  de- 
creasing volume  suggests  that  this  may  bo  tho  cause  of  tho 
negative  expansion  of  vitroous  silica  at  low  temperature,  A 
more  detailed  investigation  of  this  is  planned  in  a later  paper. 
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